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IS.  SURRLEMENTARY  NOTCS 


IS.  KEY  W6R0S  (Continue  on  tore too  el  do  II  nocoeoory  end  Identity  By  Block  numkor) 

Electromagnetic  Fields  Vertical  Electric  Dipole 

Horizontal  Electric  Dipole  Vertical  Magnetic  Dipole 

Horizontal  Magnetic  Dipole  Subsurface-to-Subsurface  Propagation 


10.  ABSTRACT  (Continue  on  roootoo  tide  II  nocoooocr  end  Identity  By  Block  numBec) 

New  formulas  for  the  electric  and  magnetic  fields  produced  by  the  four 
elementary  dipole  antennas  have  been  developed  for  the  subsurface-to-subsurface, 
subsurface-to-surface,  surface-to-subsurface,  and  surface-to-surface  propagation 
cases.  These  formulas  are  of  rather  simple,  but  useful,  form  and  are  completely 
general  (i.e.,  the  air  can  easily  be  replaced  by  the  sea  bottom).  They  are 
valid  at  any  frequency  and  at  any  range  beyond  a  certain  minimum  distance  for 
the  flat-earth  case.  The  main  restrictions  on  these  formulas  are  the  square  — 
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pof  the  index  of  refraction  is  VlO  and-fS^  the  horizontal  separation  is  \>3  times 
'  the  sum  of  the  depths  of  buriar  of  the  transmitting  and  receiving  point  sources. 
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With  these  new  formulas,  computer  evaluation  can  be  reduced  to  fractions 
of  a  minute,  compared  with  hours  for  the  complete  numerical  evaluation  of  the 
exact  Sommerfeld  integrals. 


There  also  will  be  an  interference  pattern  set  up  under  certain  close-range 
conditions  because  the  three  waves  (direct,  modified  mirror  image,  and  lateral) 
may  interfere,  either  constructively  or  destructively,  with  each  other. 
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GLOSSARY  OF  SYMBOLS 
D  (p2  +  h2)1/2  (meters) 

Ep  Horizontal  electric-field  component  in  the  p  direction  (volts/meter) 

Horizontal  electric-field  component  in  the  $  direction  (volts/meter) 

Ez  Vertical  electric-field  component  (volts/meter) 

F  F(a>0),  Sommerfeld  ground-wave  attenuation  factor 

h  depth  (h  0)  of  transmitting  antenna  with  respect  to  the  earth's 

surface  (meters) 

Hp  Horizontal  magnetic-field  component  in  the  p  direction 

(amperes/meter) 

tyj,  Horizontal  magnetic-field  component  in  the  (f>  direction 

(amperes/meter) 

Hz  Vertical  magnetic-field  component  (amperes/meter) 

HED  Horizontal  electric  dipole 

HMD  Horizontal  magnetic  dipole 

I  Current  (amperes) 

IQ  Modified  Bessel  function  of  the  first  kind  and  order  zero 

Ij  Modified  Bessel  function  of  the  first  kind  and  order  one 

Jq(ap)  Bessel  function  of  the  first  kind,  order  zero,  with  argument  Ap 
Kq  Modified  Bessel  function  of  the  second  kind  and  order  zero 

K1  Modified  Bessel  function  of  the  second  kind  and  order  one 

m  Magnetic -dipole  moment  (ampere -meters2) 

n  Yi/Y0,  index  of  refraction 

Y  Y 

N  IQ  -—[Rj  -  (z  +  h)]  Kq  —  [Rj  +  (z  +  h)]  =  Foster  integral 

p  Electric-current  moment  (amperes-meters) 

PQ  exp(-y1R0)/R0  =  Sommerfeld  integral  (meters-1) 

iii. 
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GLOSSARY  OF  SYMBOLS  (Cont'd) 


^PC-YjRjJ/Rj  =  Sommerfeld  integral  (meters-1) 
(p2  +  z2)1/2  (meters) 

[p2  +  (z  -  h)2]1/2  (meters) 

[p2  +  (z  +  h)2]1/2  (meters) 

Time  (seconds) 

16I1K1  +  r VCIjKj  -  I0K0)  +  4yp(I1K0  -  IoKl) 
(*2  +  Yq)1/2  (meters-1)  (air) 

(X2  +  y2)1/2  (meters-1)  (earth) 

Vertical  electric  dipole 
Vertical  magnetic  dipole 
31^  -  (yp/2)(I0K1  -  IjKg) 

Sommerfeld  numerical  distance 


Depth  (z  >_  0)  of  receiving  antenna  with  respect  to  earth's 
surface  (meters) 

(-a)2u0eo^ 1/2 »  upper  half- space  (free-space)  propagation 
constant  (meters-1) 


(ioHJjOj  -  (o2PjE j) 1  /2 ,  lower  half-space  (earth)  propagation 
constant  (meters'1) 


-1) 


1/2 


me. 


’1  J 


-1/2 


skin  depth  in  the  lower 


half-space  (earth)  (meters) 


10"9/36if  farads/meter,  permittivity  of  free  space 


Permittivity  of  lower  half-space  (earth)  (farads/meter) 

Dummy  integration  variable  in  the  basic  Sommerfeld 
integrals  (meters-1) 

(x2  +  y2)1/2  radial  distance  in  s  cylindrical  coordinate 
system  (meters) 

Conductivity  of  the  lower  half-space  (earth)  (Siemens/raeter) 


* 


IV 


TR  6881 


GLOSSARY  OF  SYMBOLS  (Cont'd) 

4>  tan_1(y/x),  azimuth  angle  in  a  cylindrical  coordinate  system 

p  =  pQ  =  4tt  x  10" 7  henries/meter,  permeability  of  free  space 
tan-1[(z  -  h)/p],  elevation  angle 
ipj  tan_1[(z  +  h)/p],  elevation  angle 

(o  2irf  radians/second,  angular  frequency 
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NEW  FORMULAS  FOR  HED,  HMD,  VED,  AND  VMD 
SUBSURFACE-TO-SUBSURFACE  PROPAGATION 


INTRODUCTION 


In  two  recent  papers,  Wu  and  King* >2  have  derived  new  simple  formulas  for 
the  electric-field  components  generated  by  a  horizontal  electric  dipole  (HED) 
in  a  half-space  of  water  or  earth  near  its  boundary  with-  air.  They  claim  that 
their  formulas  are  valid  for  p2  >>  (z  +  h)2  and  Jn2|  »  1,  where  n  (=  Yj/Yq)  is 
the  index  of  refraction.  On  examining  their  results,  the  author  of  this  report 
has  discovered  that  they  could  have  been  obtained  almost  by  inspection  from  the 
previously  derived  results  of  Wait, 3  Weaver,**  Bannister, 5  and  Bannister  and 
Hartb  (most  of  which  are  summarized  in  Kraichman?) . 

In  the  past,  many  investigators  erroneously  have  believed  that  the  field- 
strength  equations  tabulated  in  chapter  3  of  Kraichman  are  only  valid  when  the 
conduction  currents  in  the  water  or  earth  are  much  greater  than  the  displace¬ 
ment  currents  (i.e.,  Oj  »  ojEj).  Indeed,  as  long  as  |n2|  »  1,  the  displace¬ 
ment  currents  can  be  included  simply  by  replacing  Oj  with  +  iu^  in  the 
field- strength  equations.  Thus,  Kraichman 's  tabulated  results  are  considerably 
more  general  than  they  are  stated  to  be. 

It  is  the  purpose  of  this  report  to  present  new  formulas  for  HED.  hori¬ 
zontal  magnetic  dipole  (HMD),  vertical  electric  dipole  (VED),  and  vertical 
magnetic  dipole  (VMD)  subsurf ace- to-subsurf ace  propagation.  These  formulas 
have  been  obtained  completely  from  previously  derived  results.  The  main 
restrictions  on  their  use  are  (1)  the  square  of  the  index  of  refraction  is 
>10  and  (2)  the  horizontal  separation  is  greater  than  or  equal  to  three  times 
the  sum  of  the  depths  of  burial  of  the  transmitting  and  receiving  point  sources. 
An  additional  restriction  must  also  be  applied  to  the  lateral-wave  components. 
The  quantity  |yjP2/(z  +  h) |  must  be  greater  than  or  equal  to  4Cj,  where  cx  = 

3,  6,  9,  15,  or  25,  depending  on  the  particular  field-strength  component.  This 
restriction  also  applies  to  Wu-and-King's  results* > 2  and  to  those  tabulated  in 
tables  3.7  and  3.16  of  Kraichman. 7  These  new  formulas  also  avoid  the  use  of 
the  unphysical  distance  z  +  h  +  p  employed  by  Wu  and  King .  *  >  2 

In  this  report,  the  four  dipole  antennas  (VED,  VMD,  HED,  and  HMD)  are 
situated  at  depth  h  (h  >_  0)  with  respect  to  a  cylindrical  coordinate  system 
(p,<j>,z)  and  are  assumed  to  carry  a  constant  current,  I.  The  axes  of  the  VED 
and  HED  (of  dipole  moment  p)  are  oriented  in  the  z  and  x  directions,  respec¬ 
tively,  and  the  axes  of  the  VMD  and  HMD  (of  dipole  moment  m)  are  oriented  in 
the  z  and  y  directions,  respectively.  The  earth  or  water  occupies  the  lower 
half-space  fz  >  0)  while  the  air  occupies  the  upper  half-space  (z  <  0).  The 
magn^-ic  permeability  of  the  earth  is  assumed  to  equal  p0,  the  permeability  of 
free  space.  Meter-kilogram-second  (MKS)  units  are  employed  and  a  suppressed 
time  factor  of  exp(iwt)  is  assumed. 
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WU-AND-KING'S  METHOD 


In  their  first  article,!  Wu  and  King  developed  a  new  simple  and  very 
accurate  formula  for  the  radial  electric  field  (Ep)  of  a  HED  in  a  dissipative 
or  dielectric  half-space  near  its  boundary  with  air.  They  examined  in  detail 
the  interference  patterns  generated  by  the  direct  and  lateral  waves  that  orig¬ 
inate  at  the  dipole  for  three  values  of  £j ,  numerous  values  of  Oj,  and  a  wide 
range  of  frequencies.  They  confirmed  the  accuracy  of  the  new  formula  by  com¬ 
parison  with  numerically  evaluated  Sommerfeld  integral  results. 

In  their  second  articled  Wu  and  King  developed  new  simple  formulas  for 
the  transverse  (E^)  and  vertical  (Ez)  electric-field  components  generated  by  a 

buried  HED  source.  They  then  compared  all  three  electric-field  components 
with  numerical  integration  results  for  the  case  where  a1  =  3.5  S/ra,  ex  =  45e0, 
and  f  =  600  MHz  (see  figure  1  of  Wu  and  King2) .  The  agreement  between  the 
simple-formuia  and  numerical-integration  results  was  excellent  for  the  radial 
component  when  |y,p|  >  2.8  and  p  >  2(z  +  h).  However,  substantial  agreement 

between  the  simple-formula  and  numerical -integration  transverse  and  vertical 
components  was  not  achieved  until  j y 1 p |  >  11.5  and  p  >  8(z  +  h) .  Furthermore, 
the  simple  formulas  predicted  a  dip  in  the  transverse  component  and  no  dip  in 
the  vertical  component  near  lyjpj  ~  8,  while  the  numerical-integration  results 
predicted  the  opposite. 

Because  the  radial -component  simple  formula  predicts  the  interference 
patterns  very  accurately  and  the  transverse-  and  vertical -component  simple 
formulas  do  not,  it  is  apparent  that  Wu  and  King  have  made  some  errors  in 
their  transverse-  and  vertical -component  derivations.  A  prime  suspect  is  the 
unphysical  distance  term,  z  +  h  +  p,  which  appears  in  the  E,  and  Ez  formulas 
but  does  not  appear  in  the  Ep  formula.  Wu  and  King  could  also  have  made  a 
sign  error  so  that  the  direct  and  lateral  waves  are  adding  instead  of  sub¬ 
tracting...  or  vice  versa. 

The  procedure  employed  by  Wu  and  King! , 2  caR  best  be  shown  by  example. 
When  both  the  transmitting  and  receiving  dipoles  are  located  below  the  earth's 
surface  (h  and  z  >_  0) ,  the  HED  nx  vector  may  be  expressed  as  (Wait3) 


n  -  7 — .  ^ — : - c-[Pn  -  P,  + 

x  4ir(o1  +  i(i)E1) 1  0  1 


where 


is  the  direct -wave  contribution  and 


(2) 


(3) 


2 
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is  the  mirror-image  contribution.  The  remaining  term  is 

-uj (z+h) 


f  e-u'V 

=  21  - - Jn(Ap)AdA 

J  uj  +  u0  O'- 

0 


O' 

-  ypf0  Cu‘  ■  u°,e  1 


(4) 


U+h) 


Jn(Xp)XdA 


where 


=  p2  +  (z  -  h)2. 


R2  =  p2  +  (z+h)2. 


2  -  o  o 

un  =  X2  + 


u.  =  X2  +  y 


2 

1* 


Uj  =  -io2yoeo>  and 


yJ  =  iuy0(ai  +  iwe^. 


For  | n2 |  »  1  and  p2  »  (z+h)2,  Wu-and-King’s  procedure  is  to  set  RQ 
and  Rx  equal  to  p  everywhere  except  in  the  exponents  and  let 


-ux(z+h)  -Yj(z+h) 


-  e 


(5) 


in  equation  (4).  Therefore, 


-TlRo 


P0  ' 


P 

-YlRl 


(6) 

(7) 


and 


I  - 


2e 


-Yj (z+h) 


Ct? 


I  2  v"  /• 

2  " ~  J  -  U 
l  V  Jn 


0)J0(Ap)\dA  . 


(8) 


This  integral  can  be  readily  evaluated  (Wait,3  Erdelyi8)  to  yield 


;  3 
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2e"Y‘lZ*h'  f  -T,o 

(Y?  -  Y§)P^'  *  Y‘P)e 


-  (1  +  Y0P)e 


•Y0p 


]• 


(9) 


Thus,  by  following  the  procedure  of  Wu  and  King,l>2  the  final  expression 
for  nx  is  (for  |n2|  »  1  and  p2  »  (z  +  h)2) 


p 

L-y.Ko 

e-YA 

4ir(a1  +  iuiEj) 

p 

p 

1  ,  (10) 
2  ,  -Yi(z+h+p)  ^  %  -YnP  “Yt  (z+h)"|  [ 

•  [t2T  y2)V3  LQ  *  YlP)e  e  1  JJ  • 

Note  that  IIX  has  four  components:  (1)  a  direct-wave  component,  (2)  a  mir¬ 
ror-image  component,  (3)  a  lateral -wave  component,  and  (4)  a  false  component 
that  depends  on  the  unphysical  distance  z  +  h  +  p. 


BANNISTER'S  METHOD 


As  we  shall  soon  see,  the  unphysical  distance  z  +  h  +  p  can  be  avoided. 
When  the  measurement  distance  is  much  less  than  a  free-space  wavelength,  equa¬ 
tion  (1)  reduces  to 


Pn  —  P_  +  _ 

0  1  (t2 


H-r/  c», 
i  -  rPl  ' 


>Je 


-Uj (z+h) 


x  4-rrCaj  +  iuCj) 

Wait3  has  shown  that  this  equation  is  equivalent  to 

nv  - 


J  0 (Xp)XdX 


p 

P  P  i  ^ 

(d2l>1  .  33N  +  y2  3n\ 

4n(a1  +  iuEj) 

[°  1  (Y?  -  Yq) 

^  9z2  3z^  1  3z  y 

where 


N  =  I, 


2  (Ri  "  Cz  +  h)] 


K0  2  [R1  +  (Z  +  h)] 


Wait3  has  also  shown  that,  when  |yxp 1  »  1  and  p  »  (z  +  h), 

-Yj (z+h) 


N  - 


YiP 


(ID 


(12) 


(13) 


(14) 


Quoting  Wait,3  "Hie  manner  in  which  the  exponential  factor  exp[-Yj(z  +  h)] 
occurs  is  rather  interesting.  It  is  only  in  the  integral  N  that  this  factor 
emerges . " 
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Employing  equation  (14)  and  taking  the  indicated  derivatives  in  equation 
(12)  results  in 


•Yi(z+h)  -YjRj 


(y2  -  Yq)p3  ri 


p  (e~YlR0  ,  2e~Yl(Z+h) 

4tt(o1  +  itue^  j  R0  +  (y2  _  Y2)p3  '  ^ 

”^1R1  1 

(yf-  rg)»f  f(1  '  3  ^  *  TlRl)  "  Y*R‘  Si”2  *‘4' 


2e-nRi 


where  sin  if>  =  (z  +  h)/R1. 

For  p2  »  (z  +  h)2,  we  can  set  RQ  and  R1  equal  to  p  everywhere  except  in 
the  exponents.  Thus,  1 

-Y,R0  _  -Y,(z+h) 

n  -  -  e-1-0  ,  2e  - . 

x  4ir(o1  +  iuEj)  |  p  (y2  -  y2)p3 

-  +  (1  *  YlP)j  J  • 


(yf  -  Yq)p2 


Here,  we  see  that  nx  has  only  three  components:  (1)  a  direct-wave  component, 
(2)  a  modified  mirror-image  component,  and  (3)  a  lateral-wave  component. 

There  is  no  false  component  that  depends  on  the  unphysical  distance  z  +  h  +  p. 

At  first  glance,  this  procedure  appears  to  be  considerably  more  compli¬ 
cated  than  the  procedure  employed  by  Wu  and  King.  Luckily,  however,  the  gen¬ 
eral  quasi-static  range  field  components  have  already  been  derived  for  the 
four  elementary  dipoles  (Wait, 3  Weaver, 4  Bannister, 5  Bannister  and  Hart. 6  Wait 
and  Campbell, 9.10  and  Sinya  and  Bhattacharyall)  as  have  the  quasi-near,  near¬ 
field,  and  farfield  range  lateral -wave  expressions  (tables  3.2,  3.7,  and  3.16 
of  Xraichman7) .  Although  displacement  currents  were  ignored  in  most  of  these 
analyses,  they  can  be  included  simply  by  replacing  Oj  by  0j  +  ime^  (as  long  as 

|n2|  »  1).  For  convenience  sake,  the  direct  and  modified  mirror- image  con¬ 
tribution  to  each  dipole  field-strength  component  is  listed  in  the  appendix. 

VF 

For  some  components,  these  expressions  are  of  very  simple  form  (H^  ,  equation 

(A-15)),  while  for  other  components,  these  expressions  are  quite  complicated 

ChY*4,  equation  (A-18)). 

2 

Thus,  the  hardest  part  of  the  problem  has  already  been  solved.  We  can 
now  use  these  previously  derived  results  to  obtain  adequate  formulas  for  the 
fields  produced  by  submerged  dipole  sources  subject  to  the  conditions  |n2[  » 

1  and  p2  »  (z  +  h)2. 

The  derivation  procedure  that  we  will  follow  is  (1)  take  the  previously 
derived  direct  and  modified  mirror-image  results  (see  appendix)  and  let  p  be 
>_  3(z  +  h),  remembering  not  to  replace  R0  and  Rj  by  p  in  the  exponents,  and 

(2)  for  the  lateral -wave  expressions,  let 


f(p,z  +  h)  =  f(p,0)e 


-Yj(z+h) 
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where  0  refers  to  an  infinitesimal  distance  below  the  surface  of  the  earth  or 
water. 

The  quasi-static  range  (|y0p|  <<  l)  lateral-wave  functions  ft.p,0)  can  be 
obtained  from  Bannister^  or  table  3.8  of  Kraichman?  (with  o1  replaced  by  Oj  + 
iucj) .  For  convenience  sake,  the  lateral-wave  functions  f(p,0)  exp[-Yj(z  +  h)] 
are  presented  in  table  1.*  Note  that,  for  some  of  the  field-strength  compo¬ 
nents,  these  expressions  are  identical  to  the  quasi-near  range  (|y0p|  <<:  1 
«  | yxP 1 )  results  (see  table  3.16  of  Kraichman7) .  Also,  note  that  some  of  the 
signs  are  different  than  those  of  Bannister^  and  Kraichman.7  This  is  because, 
in  this  report,  we  have  inverted  the  coordinate  system  so  that  z  and  h  are 
positive  depths. 

Approximately  half  of  the  field-component  formulas  listed  in  table  1 
involve  products  of  modified  Bessel  functions  of  argument  YiP/2.  Numerical 

values  for  these  functions  have  been  provided  by  Bannister. 5  when  Iy^I  £.4, 
the  function  Y^pl^  *  1,  while  when  | y jP  I  >  6,  the  function  yxpw  -  2.  Fur¬ 
thermore,  when  | YiP |  1  10,  the  function  YiPT/2  ~  3. 

As  an  example  of  our  derivation  procedure,  consider  the  HED  radial  elec¬ 
tric  field-strength  component.  The  quasi-static  range  lateral-wave  Ep  compo¬ 
nent  can  be  obtained  from  table  1  or  from  equation  (67)  of  Wait. 3  it  is  equal 
to 


ehe _ p  cos  $  e~y  i  tz+h) 

p  2ir(0j  +  iwe^p3 


(18) 


The  nearfield  range  lateral-wave  Ep  component  can  be  obtained  from  table 
3.7  of  Kraichman?  (with  Oj  replaced  by  o1  +  iioCj),  or  from  equations  (23)  and 
(32)  of  Wait. 3  Therefore, 


-Y,(z+h)  -y0P 

eHE  _  p  cos  <fre  *  e  u 
p  2ir(o1  +  iuie^p3 


(1 


+  YnP2) 


(19) 


The  farfield  range  lateral-wave  Ep  component  can  be  obtained  from  table 
3.2  of  Kraichman?  (with  Oj  replaced  by  Oj  +  iuiCj),  or  from  equations  (13)  and 
(32)  of  Wait. 3  Thus, 


eHE 

“p 


-Y,(z+h) 

p  cos  4>e  _ e 

2ir(a1  +  iioe^p3 


"V 


(YqP2F)  » 


(20) 


where 


F  =  F(wq)  =  1  -  i(irw0)1/2e  0erfc(iw01/2)  (21) 

is  the  Sonmerfeld  surface-wave  attenuation  function  and 

‘Tables  have  been  placed  together  at  the  end  of  this  report. 
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w 


0 


V 


2n2 


(22) 


is  the  Sommerfeld  numerical  distance.  For  small  numerical  distances  F(wQ)  -  1, 
while  for  large  numerical  distances  F(w0)  -  -l/(2w0). 


When  |n2|  »  1,  the  range  of  validity  of  equations  (18)  and  (19)  overlap 
when  |y0p|  «  1.  Similarly,  when  jwQj  «  1  and  (y0p|  »  1,  the  range  of  valid¬ 
ity  of  equations  (19)  and  (20)  overlap.  Therefore,  we  can  simply  combine  equa¬ 
tions  (18) ,  (19) ,  and  (20)  to  obtain  an  expression  for  the  Ep  lateral -wave 
component  valid  from  the  quasi-static  to  the  farfield  range.  Therefore,  for 
|n2|  »  1  and  p2  »  (z  +  h)2, 


eHE  _  P  cos  4>e 
p  2tt(o 


-Yj(z+h)  -y0p 


imepp' 


(1  +  YnP  +  YnP2p) 


(23) 


The  direct  and  modified  mirror-image  contribution  to  each  dipole  field- 
strength  component  is  listed  in  the  appendix.  These  expressions  are  valid  for 
|n2[  •»  1.  When  p  >  3(z  +  h),  equation  (A-l)  reduces  to 


PHE 

CP 


p  cos  4> 

2ir(o1  +  iwe^p3 
(3  +  3YlP  +  Y2P2 


(1  +  Y2P)e  YlR° 


(Z  +  hj^-Y^! 


2p' 


(24) 


We  can  now  combine  equations  (23)  and  (24)  and  obtain  an  expression  for 
the  HED  Ep  component  valid  at  almbst  any  range  from  the  source  subject  to  the 
conditions  that  |n2|  »  1  and  p  ^  3(z  +  h) .  The  final  expression  is 


pHE  _ 
EP 


p  cos 


(1  +  Y0P  +  YoP2F)e  Y°Pe  Yl("Z  +  (1  +  YxP)e  YlR° 


2tt(o1  +  irne^p3 

rt  .  t  .  ?  p-iT  (z  -  h)2  (z  +  h)2  "Y]®]]  ( 

-  (3  +  3YlP  +  Y2P2)[A-— 2'^  o  +  ■■  2p2~'~~c  j)  * 


(25) 


When  p2  »  (z  +  h)2,  the  last  two  terms  of  equation  (25)  a.<*  negligible 
compared  with  the  first  two,  resulting  in 


CHE 


p  cos  *  L  9  9^  "TnP  *Y,  (z+h 

Tt K1  +  T0P  +  YpPzF)e  0  e  1 
+  ioie  )p3[  0  0 


+  (1  +  YjP)e 


‘y1R0 


(26) 


which  is  identical  to  Wu-and-King's  result *  for  |n2i  »  1  and  p2  »  (z  +  h)2. 

As  we  mentioned  previously,  Wu  and  King  have  shown  that  this  simple  formula  is 
in  excellent  agreement  with  the  exact  Sommerfeld  integral  numerical-integration 
results. 
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New  formulas  for  the  electric  and  magnetic  fields  produced  by  the  four 
elementary  dipole  antennas  are  presented  in  tables  2  through  9  for  the  subsur- 
face-to-subsurface,  subsurf ace-to-surface,  surf ace-to-subsurf ace,  and  surface- 
to-surface  propagation  cases.  All  of  these  formulas  have  been  obtained  from 
previously  derived  results  according  to  the  procedure  outlined  above  and  are 
strictly  valid  for  jn2|  »  1  and  p2  >>  (z  +  h)2.  (However,  for  most  cases, 
the  requirement  that  |n2|  10  and  p  ^  3(z  +  h)  is  sufficient.) 

It  should  be  noted  that  for  many  (but  not  all)  cases,  the  range  of  valid¬ 
ity  of  the  formulas  presented  in  tables  2  through  7  can  be  extended  down  to 
p  -  (z  +  h)  if  the  direct  and  modified  mirror-image  terms  in  these  equations 
are  replaced  with  the  equations  listed  in  the  appendix.  For  example,  the  HED 
Ep  expression  (equation  (25))  would  be  replaced  by 


EHE  ,  _ E 

P  4tt(cj 


cos  4>  2  2  2ci  Yop  Y^fz+h) 

- ■  — ^ — s-(l  +  Y*P  +  YnP^F)®  u  e  1 

x  +  loiCj)  p3v  '(r  '<r  1 


+  [(3  cos2  <f>0  -  1)(1  +  YjRq)  -  y2R2  sin2  ^0]— 


'l"0 


-  (3  +  3y1R1  +  YjR^)  sin2  <|»1- 


where  sin  =  (z  -  h)/RQ  and  cos  ^  =  p/Rfl. 

For  the  subsurface-to-surface  and  surface-to-surface  propagation  cases 
(tables  4  and  8),  the  vertical  electric-field  (Ez)  receiving  antenna  is 
assumed  to  be  located  an  infinitesimal  distance  above  the  earth's  surface.  To 
obtain  expressions  for  the  vertical  electric  fields  just  below  the  surface, 
multiply  the  Ez  equations  in  tables  4  and  8  by  1/n2. 


COMPARISON  OF  BANNISTER  AND  WU-AND-KING  RESULTS 


In  attempting  to  explain  the  major  discrepancy  between  the  HED  E^  and  Ez 
simple-formula  and  numerical -integration  results  near  |yjP |  ~  8  of  their  fig¬ 
ure  1,  Wu  and  King2  stated  that  very  small  changes  in  frequency  significantly 
alter  the  interference  pattern  so  that  close  agreement  in  a  small  range  near 
such  a  region  cannot  be  expected.  This  statement  is  in  direct  opposition  to 
the  Ep  results  presented  in  their  first  article  (Wu  and  King1) ,  where  sub¬ 
stantial  agreement  between  the  simple-formula  and  numerical -integration  results 
was  achieved  in  the  interference  region  when  |n2|  »  1  and  p  ^  5(z  +  h). 

They  also  noted  that,  at  greater  distances  where  the  lateral  wave  domi¬ 
nates,  the  Ep,  E^,  and  Ez  expressions  were  highly  accurate.  They  certainly 
should  be  highly  accurate  at  the  greater  distances  because  their  lateral-wave 
formulas  are  essentially  equivalent  to  Wait's  results, 3  which  have  been  suc¬ 
cessfully  utilized  for  over  20  years. 


&4SLSi2Ss>Ha2£ 


'v£*\*«lfWTOrr<eM'-' 


ran? 
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The  major  value  of  Wu-and-King's  work  is  in  determining  adequate  field- 
strength  expressions  in  the  range  where  the  lateral  and  direct  waves  interfere 
(either  constructively  or  destructively) .  They  have  succeeded  with  the  HED  Ep 
component.  However,  they  have  made  some  errors  in  their  derivation  of  the  HED 
E^  and  Ez  components. 

The  HED  E^  and  Ez  expressions  derived  in  this  report  (table  2)  are  in  very 
good  agreement  with  the  numerical -integration  results  presented  in  figure  1  of 
Wu  and  King2  when  p  >_  3(z  +  h) .  They  correctly  predict  a  substantial  dip  in 
the  E7  component  (due  tc  interference  between  the  direct,  mirror-image,  and 
lateral  waves)  near  |  ~  8  and  no  dip  in  the  E.  component  at  this  range 
(which  is  directly  opposite  to  Wu-and-King's  simple-formula  results). 

To  see  where  Wu  and  King  erred,  we  will  compare  their  E^  and  Ez  formulas 

with  the  expressions  listed  in  table  2  for  the  situation  where  the  Sommerfeld 
numerical  distance  is  small  (i.e.,  F  -  1) ,  p  ^  S(z  +  h),  and  | yxP |  >,4.  For 
this  situation,  the  E^  and  Ez  expressions  in  table  2  reduce  to 

EHE  -  ?.»*.♦■  L  .  -  (1  * 

v  2ir(o1  +  lme^p3  [ 

/  v 1  (28) 


HE - fa>?t  ,  p)e'YoPe‘Yltz*h) 

2tt  (a ~  +  iuje^p2  n2  0 

-  ~^j(3  +  3YjP  +  yfp2)  [(z  -  h)e  YlR°  +  (z  +  h)e  YlR! 


When  F~l,  p>_S(z+h),  and  j  y  Lp  (  >_  4,  Wu-and-King's2  E^  and  Ez  expres¬ 
sions  reduce  to 

EHE  .  P  sin  <j>  2f,  D-,  "V  "YlCz+h) 

E*  2tt(Oj  7  ime^psp1  W  6 

-  (2  +  2yjP  +  Y2P2)e  Yl*"  +  (1  +  YjP  <•  YjP2)©  YlRl  (30) 

+  j(l  +  YXP  +  Y2P2)^e  YlR°  -  e  YlRl^j 


EHE  .  -^L-SP-LA-  !i(1  +  Y0p)e'Yope-Yi(z+h) 

2ti(o1  +  ime^p2  n2  0 

+  -^y(3  +  3yxp  +  y2p2)  [(z  -  h)e  YlR°  +  (z  +  h)e  YjRl 
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A  comparison  of  the  two  expressions  (equations  (28)  and  (30))  reveals 
that  the  difference  between  them  is  the  unphysical  distance  (z  +  h  +  p)  term 
of  Wu  and  King.  If  we  replace  z  +  h  +  p  by  the  physical  distance  Rlf  Wu-and- 
King's  formula  reduces  to 


f  -  ,  --3  2d  * 

9  2ir(o1  +  i<oe1)pd  [  u 

+  +  Yjp  +  yfp^^e  YlR°  -  e  YlRlj  , 


(32) 


which  is  identical  to  equation  (28) . 


A  comparison  of  the  two  Ez  expressions  (equations  (29)  and  (31))  reveals 
that  the  reason  that  Wu-and-King ' s  formula  does  not  predict  the  dip  in  the 
interference  region  is  because  of  a  sign  error.  Their  results  indicate  that 
the  sum  of  the  direct  and  mirror-image  waves  add  to  the  lateral  wave  whereas, 
in  reality,  the  sum  subtracts  from  the  lateral  wave. 


RANGE  OF  VALIDITY  OF  LATERAL-WAVE  FORMULAS 


As  we  have  previously  mentioned,  the  quasi-static  range  lateral-wave 
exponential  attenuation-with-depth  factor  exp [-’Yj  (z.  +  h)]  emerges  only  from 
the  integral  N,  where 


N  = 


"2~[R1  -  (z  +  h)j 


+  (z+h)] 


(13) 


When  |YjP|  »  1  and  p  »  (z  +  h),  Wait  and  Campbell^  have  shown  that  the 
modified  Bessel  functions  may  be  replaced  by  only  the  first  terms  in  their 
respective  asymptotic  expansions  to  obtain  adequate  HMD  quasi-near  range  hori¬ 
zontal  magnetic-field-component  expressions.  On  the  other  hand,  Sinha  and 
Bhattacharyall  have  shown,  for  the  VMD  case,  the  first  two  terms  of  the  modi¬ 
fied  Bessel  function’s  asymptotic  expansions  must  be  employed.  This  indicates 
that  the  range  of  validity  of  the  quasi-static  range  lateral-wave  expressions 
will  not  be  the  same  for  all  field  components. 


As  an  example,  consider  the  quasi-static  range  HED  Ep  lateral-wave  compo¬ 
nent,  which  can  be  expressed  as 


p  cos  j> 


2-n(.Oi  +  iuu  ,)p  3P 


00 

l 


-Uj (z+h) 

Jr 


(Xp)dX 


(33) 


When  IyjpI  »  1  and  p  »  (z+h),  the  usual  procedure  is  to  replace  in  the 
exact  integral  expressions  by  y j  >  the  propagation  const.?  <t  in  the  earth, 
resulting  in 
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cLW  -  P  cos 


-YjCz+h)  „  r~ 

- t — Y-l  J.  (Xp)dX 

j  +  luje^P  3p  J  O'- 


Since  this  integral  is  equal  to  1/p  (Erdelyi8) , 


,-Yx(z+h)  -Y.(z+h) 

- - -^=f(p,0)e,i 


-LW  _  P  cos  <fre 


p  217(0^  ime-j^p3 

So  far,  we  do  not  know  exactly  what  values  |yxP|  and  p/ (z  +  h)  must  have 
in  order  to  utilize  equation  (35)  or  any  of  the  other  quasi-static  range  lat¬ 
eral-wave  formulas  presented  in  table  1. 


As  a  first  order  approximation,  we  will  let 
ux  =  (X2  +  y| ) 1/2  ~  Yi  + 


so  that 


-Uj(z+h)  _  e-Y1(z+h)e-X2(z+h)/2Y1  „  ^YjU+h)^  _  X2(z  +  h)  j 


Inserting  equation  (37)  into  equation  (33)  results  in 


-Y.Cz+h) 

]  cos  j>e  1 _ £_ 

2tt(Oj  +  io>ei)p  3p 


J0(Xp)dX 


o 

(z  +  h)  r 

2Y1  l 


X2Jn(Xp)dX  . 


Since  the  second  integral  is  equal  to  -1/p3  (Erdelyi8), 


-YjCz+h) 

eLW  „  P  cos  $6 

p  2ir(o1  +  itoe1)p3 


3  1  +  2y 


J_hl 

1’2  I 


It  can  easily  be  shown  that  the  error  incurred  in  neglecting  the  second 
term  is  less  than  1  dB  if  the  quantity 


frlP^z  +  h^  - 


4Cj  , 


where  cy  =  3  for  the  HED  Ep  component.  That  is,  to  a  first  order  approximation, 
it  is  the  product  of  Iy^pI  times  p/ (z  +  h)  that  must  exceed  a  specified  number 
to  accurately  utilize  equation  (35)  or  any  of  the  other  quasi-static  range 
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lateral-wave  formulas  presented  in  table  1.  For  example,  if  p  =  5(z  +  h) ,  the 
quantity  IyjpI  must  be  >2 .4  to  use  equation  (35). 

The  most  severe  restriction  will  be  for  the  VMD  Hz  component.  Following 
the  same  procedure  as  in  the  derivation  of  the  HED  Ep  component,  we  can  express 
the  VMD  Hz  component  as 


-Y  (z+h) 
9me  1 


2tt  CYf 


yJ)p5 


25  (z  +  h] 
2y,P2 


Again,  it  can  easily  be  shown  that  the  error  incurred  in  neglecting  the  second 
term  is  less  than  1  dB  if  the  quantity 


where  Cj  =  25  for  the  VMD  Hz  component.  For  example,  if  p  =  5(z  +  h),  the 
quantity  | y XP i  must  be  >20  to  employ  the  quasi-static  range  VMD  Hz  lateral-wave 
formula  presented  in  table  1. 


The  values  of  Cj  and  range  where  f(p,z  +  h)  can  be  replaced  by  f(p,0) 
expt-YjCz  +  h)]  (i.e.,  the  range  where  the  quasi-static  range  lateral-wave 
formulas  presented  in  table  1  can  be  used)  are  presented  in  table  10  for  each 
field-strength  component.  Here,  we  see  that  for  the  HED  and  HMD  Ep,  Ez,  and 
H^  components  and  the  VED  Ep  and  H^  components  (cl  =  3),  the  quantity  1  y XP2/ 

(z  +  h) |  must  be  >12,  while  for  the  HED  Hz  and  VMD  E^  and  Hp  components  (cx  = 
15),  the  quantity  |y1P2/(z  +  h) |  must  be  >60. 

This  restriction  (equation  (40))  also  applies  to  Wu-and-King's  results1'2 
and  to  the  quasi-near  and  nearfield  range  subsurface-to-subsurface  propagation 
equations  tabulated  in  Kraichman.? 

To  a  first  order  approximation,  the  range  of  validity  of  the  equations 
listed  in  table  1  can  be  extended  by  multiplying  the  field-component  expres¬ 
sions  by  the  quantity 


1  + 


Cj(z  +  h) 

“IT?”-’ 


(42) 


where  the  value  of  Cj  for  each  component  is  given  in  table  10.  For  example, 
the  HED  Hz  component  listed  in  table  1  is 


,HE  3p  sin  <fe 


-YxCz+h) 


2if(y\  ~  Y2)p4 


Since  Cj  =  15  for  this  component  (table  10),  then 
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H! 


,HE  ,  3p  sin  «j»e 


-Yj(2+h) 


2*(y2  -  Yq)p1 


x  +  15(z  Oi) 
l<  “ 


2YiP2 


(44) 


ANALYTICAL  CONFIRMATION  OF  FRASER-SMITH 
AND  BUBENIK  VMD  Hz  NULL 


A  few  years  ago,  Fraser-Smith  and  Bubenikl2  numerically  evaluated  the 
exact  Sommer f eld  integrals  and  found  a  rather  deep  null  in  the  VMD  Hz  component 
for  the  subsurface-to-surface  propagation  case  (see  figure  6  of  Fraser-Smith 
and  Bubenikl2) .  For  their  particular  situation,  the  frequency  was  100  Hz,  the 
VMD  source  depth  of  burial  was  100  m,  and  the  null  occurred  around  p  ~  250  m. 
Since  the  skin  depth,  6,  in  sea  water  at  100  Hz  is  ~25  m,  p/6  ~  10  and  h/6  ~  4. 
To  our  knowledge,  this  null  has  not  been  analytically  confirmed.  The  subsur- 
face-to-surface  VMD  Hz  component  is  equal  to  (from  table  5) 


m 


2*(Yl  -  Yg)p‘ 


(9  +  9ynp 


+  4y?P2 


+ 


(9  +  9yxP  +  4y2P2  +  YjP3) 


h2, 


- j(90  +  90YjP  +  39y2P2  +  9y3P3  +  Y4P4) 


(45) 


where  D2  =  p2  +  h2. 


For  f  =  100  Hz  and  p  =  250  m,  |y0p|  <<c  1  and  the  Y0P  terms  in  the  lateral- 
wave  portion  of  equation  (45)  are  negligible.  Furthermore,  since  the  quantity 
j YjP2/h j  ~  35,  which  is  <JC0  (see  table  10),  equation  (42)  must  be  employed. 

The  result ' ng  expression  for  Hz  is 


jVM 


-Y,h 

9me  1 
2ttY2P5  2  * 


(46) 


where 


/l  +  25h  \ 

-Y,(D-h)  ‘ 

e  1 

(  2YiP2/ 

9 

(9  +  9yjP  +  4y2P2  +  Y3P3) 


h2 

-  ~(90  +  90yxP  +  39y2P2  +  9y3P3  +  Y4P4) 


(47) 


For  p/6  -  10,  the  dominant  terms  will  be  the  y3P3  and  y4P4  terms.  To  a 
first  order  approximation,  the  y^P3  terms  will  cancel  and,  since  y4P4  = 

~4 (p/6) 4, 
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|f)2(|)'e'(D-h,/4e-i(D-h,/i  •  («) 

This  equation  will  be  at  a  minimum  near  (D  -  h)/6  =  2ir,  which  corresponds 
to  p  -  240  m.  Tlie  normalized  VMD  vertical  magnetic  field  (Hz  from  equation 
(47)  )  is  plotted  in  figure  1  versus  the  horizontal  distance,  p.  From  this  fig' 
ure,  we  can  see  that  a  rather  deep  null  (~20  dB  drop  in  field  strength  compared 
to  the  asymptotic  value)  occurs  at  a  range  of  -240  m.  This  null  is  clearly  due 
to  the  destructive  interference  between  the  direct  and  lateral  waves. 


CONCLUSIONS 


New  formulas  for  the  electric  and  magnetic  fields  produced  by  the  four 
elementary  dipole  antennas  have  been  developed  for  the  subsurface-to-subsurface 
propagation  case.  These  formulas  have  been  obtained  completely  from  previously 
derived  results.  The  main  restrictions  in  their  use  are  (1)  the  square  of  the 
index  of  refraction  is  >10  and  (2)  the  horizontal  separation  i?  >5  times  the 
sum  of  the  depth  of  burial  of  the  transmitting  and  receiving  point  sources.  An 
additional  restriction  must  also  be  applied  to  the  lateral-wave-  components. 

The  quantity  | y iP2/ Cz  +  h) j  must  be  >4c.,  where  c1  =  3,  6,  9,  15,  or  25. 
depending  on  the  particular  field-strength  component.  This  restriction  also 
applies  to  Wu-and-King ' s  recently  derived  results**2  and  to  the  subsurface-to- 
subsuvface  propagation  equations  tabulated  in  Kraichman.7 

The  range  of  validity  of  the  subsurface-to-subsurface,  subsurface-to-sur- 
face,  and  surface-to-subsurface  equations  tabulated  in  this  report  can  be 
extended  down  to  p  -  (z  +  h)  for  many  cases  if  the  direct  and  modified  mirror- 
image  terms  in  these  equations  are  replaced  with  the  equations  listed  in  the 
appendix.  The  extension  of  these  results  to  even  closer  ranges  will  be  the 
subject  of  a  future  report. 

It  should  be  noted  that  the  two  media  can  be  inverted  and  the-  air  replaced 
by  the  earth's  crust  (of  conductivity  a2  and  dielectric  constant  t2) .  The  same 
equations  (tables  1  through  10)  can  be  utilized,  as  long  as  |ii||  =  |yj/y||  i 
10  and  p  >  3(z  +  h) ,  simply  by  replacing  iwe0  by  a2  +  ime2. 

An  analytical  confirmation  of  the  Fraser-Smith  and  Bubenik  VMD  Hz  null  has. 
also  been  accoi'nlished.  This  null  is  clearly  due  to  the  destructive  inter¬ 
ference  between  the  direct  and  lateral  waves. 


-v. (z+h) 

f(p,0)e  1  Lateral  Wave  Formulas  When 
p|  «  1  [|n2|  ^.10,  p2  »  (z  +  h)2]* 
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Table  2.  Electric-Field  Subsurface-to-Sub 


Dipole 

Type 
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♦Argument  of  modified  Bessel  function  is  yxp/2. 
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irface-to-Subsurface  Propagation  Formulas  [|n2|  ^10,  p  >_  3(z  +  h)]1 
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Table  3.  Magnetic-Field  Subsurface 

1 1 n2 1  >_  10,  p 
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•Argument  of  modified  Bessel  functions  is  YjP/2. 


|Fid  Subsurface-to-Subsurface  Propagation  Formulas 
|[|n2|  >  10,  p  >  3(2  ♦  h)]* 
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Table  3.  (Cont'd)  Magnetic-Field  S 
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♦Argument  of  modified  Bessel  functions  is  YiP/2 
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♦Argument  of  modified  Bessel  functions  is  YjP/2. 
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Jo-Surface  Propagation  Formulas  [|n2|  >  10,  p  >  3h,  D  =  (p2  +  h2)1/2]* 
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Table  7.  Magnetic-Field  Surface-to-Subsurfa 
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♦Argument  of  modified  Bessel  functions  is  Yjp/2. 


Table  8.  Electric-Field  Surface-to-Surface  Propagati 


mSm 

EP 

H 

VED 

Y1P  “YnP 

— - — - -yCy  pi  K  +  Y0pF)e 

27T  (a  j  +  ill  o 

0 

VMD 

0 

ia)ji0m  r 

=  -  ,  2  ,,  u  (3  +  3yqP  + 

2ir (y2  -  YjpP  1  0 

-  (3  +  3YjP  +  Y2P2)e  YlP] 

HED 

— ~T 2 - ■  ^  •  3  [(i  +  Y0P  +  YQP2F)e  0 

2tt (cJj  +  lwe^p3  L  0  0 

+  (1  +  Y2P)e  YlP] 

7-,  p  fn  ♦  3  b  -  V<*  *  F 

27r(a1  +  lue^p*  l  0 

-  (1  +  YjP)e  YlP| 

HMD 

iup0m  cos  $ 

"  2*yiP3  (YlPllKl  +  Y0P  +  Y0p2F) 

-Y0p 

x  e  u 

iuy  m  sin  <J> 

- — [y  pw  +  y0p(i  +  f; 

21TYJP3 

♦Argument  of  modified  Bessel  functions  is  Yjp/2. 


TR  6881 


ctric-Field  Surface-to-Surface  Propagation  Formulas  (|n2|  10)* 
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Field  Surface-to-Surface  Propagation  Formulas  ( | n2  j  >  10)’ 
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Appendix 


DIRECT  AND  MODIFIED  MIRROR-IMAGE  CONTRIBUTION  TO  EACH  DIPOLE 
FIELD-STRENGTH  COMPONENT  (|n2|  »  1) 
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